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ABSTRACT 

As a part of our ongoing Volume-limited A-Star (VAST) adaptive optics survey, we have ob- 
tained observations of 26 binary systems with projected separations <100 AU, 13 of which 
have sufficient historical measurements to allow for refinement of their orbital elements. For 
each system with an estimated orbit, the dynamical system mass obtained was compared with 
the system mass estimated from mass-magnitude relations. Discrepancies between the dynam- 
ical and theoretical system mass can be explained by the presence of a previously unresolved 
spectroscopic component, or by a non-solar metallicity of the system. Using this approach to 
infer the presence of additional companions, a lower limit to the fraction of binaries, triples, 
and quadruples can be estimated as 39, 46, and 15 per cent, for systems with at least one 
companion within 100 AU. The fraction of multiple systems with three or more components 
shows a relative increase compared to the fraction for Solar-type primaries resolved in previ- 
ous volume-limited surveys. The observations have also revealed a pair of potentially young 
(< 100 Myr) M-dwarf companions, which would make an ideal benchmark for the theoretical 
models during the pre-Main Sequence contraction phase for M-dwarfs. In addition to those 
systems with orbit fits, we report 13 systems for which further orbital monitoring observations 
are required, 1 1 of which are newly resolved as a part of the VAST survey. 

Key words: binaries: general - binaries: close - techniques: high angular resolution 



1 INTRODUCTION 

The orbits of binary stars offer one of the few techniques to deter- 
mine stellar m asses, or masses an d radii in cases with a favorable 
geometry (e.g. I Andersen et aljl99lh . Orbits of binaries with young 



* E-mail: derosa@astro.ex.ac.uk 

f Based on observations obtained at the Canada-France-Hawaii Tele- 
scope (CFHT) under programme IDs 2001AF11, 2008AC22, 2009BC06, 
2010AC14, 201 1AC1 1. Based on observations made with ESO Telescopes 
at the La Silla Paranal Observatory under programme IDs 272.D-5068(A), 
074.D-0 180(A), 076.C-0270(A), 076.D-0108(A), 077.D-0147(A), 080.D- 
0348(A). Based on observations obtained at the Gemini Observatory under 
programme IDs GN2008A-Q-74, GN-2008B-Q-1 19, GN-2010A-Q-75. 



ages are particularly important to test theoretical evolutionary mod- 
els, and examples inclu de a double-lined eclipsing binary in Orion 
(e.g lStassun et al.l l2008h . Low mass stars and brown dwarfs also 
represent a regime requiring empirical calibration, and the visual 
orbits of nearby M-, L-, and T- dwarfs have been used to measure 
syste m masses and com pare with theoretical mass-luminosity rela- 
tions JPupuv et ai]|2009i) . 

Visual orbits also provide a method to search for evidence of 
additional components and determine higher order multiplicity, by 
identifying systems with dynamical masses significantly in excess 
of the theoretical predictions. These visual binaries with an indica- 
tion of unresolved companions can be monitored with spectroscopy 
or interferometry to determine the properties such as period and 
mass ratio of the closer pair and augment the statistics compiled 
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from catalogues (Tokovinin 2008). The properties of higher order 
multiple systems represent tests of formation scenario s including 
frag mentation of cores (e.g. Pringle 1989; Bonnell 2001) and disks 
(e.g. lStamatellos et al. 2007; Kra tter et al.| 2010) and may be influ- 
enced by processes such as accretion ( e.g. Bate 2000) and dynami- 
cal interactions (e.g. lMcDonald & Clarke! 1995l ; lLodato et al.l2007h 
In this paper, a subset of the systems resolved by our ongoing 
volume-limited A-star (VAST) survey are used to determine dy- 
namical system masses from orbit fits and to compare the results 
with theoretical models and search for additional unresolved stel- 
lar companions. The sample of AO-imaged binaries considered in 
this study is detailed in Section 2, and a short summary of the new 
observations is given in Section 3. The data analysis, including the 
AO image processing to determine the relative positions and the 
subsequent orbit determination from the compilation of position 
data, is explained in Section 4. Section 5 reports the astrometric 
results from the new measurements and the orbital elements and 
masses based on the orbit fits. The discussion in Section 6 covers 
a comparison with theoretical mass-magnitude relations, an assess- 
ment of the higher order multiplicity, and the identification of a set 
of targets for continued monitoring. Finally, Section 7 provides a 
summary and future directions for the project. 



2 SAMPLE 

The sample of binar ies is drawn from the ongoing VAST survey 
(De Rosa et 1), an adaptive optics (AO) imaging survey of 

A-stars within 75 parsecs, and includes the 26 systems with pro- 
jected separations less than 100 AU. The angular separations of 
the binaries range from / .'094 to 4'.' 66, and 1 1 are newly resolved. 
Figure |2]plots the measured magnitude difference as a function of 
separation for the complete sample. Of the 26 systems, 13 have a 
substantial number of previous measurements, and these systems 
comprise the orbit subsample. For the remaining 13 systems, there 
is insufficient coverage to fit an orbit, and these binaries comprise 
the monitoring subsample. Table [TJ lists each observed binary in 
the two s ubsamples, along w ith basic parameters for each such as 
distance dvan Leeuwenl|2007l) , Tycho2 Vx-band and 2MASS K a - 
band photometry teSAll 19971; Iskrutskie et alj|2006l) . and spectral 
type listed within the SIMBAD database. 

For the brightest stars within the sample, the shortest 2MASS 
exposures saturate, requiring a different method for measuring the 
photometry, resulting in significantly large r uncertainties on the 
estimated magnitudes (Skrutskie et al. 2006). Therefore for three 
of these brighter targets, near-infrared photometry was obtaine d 
from alternative sources jDucatill2002l ; iMorel & Magnenatll 19781) . 
and converted into the 2M ASS photometric bands using empiri- 
cal colour transformations jCarpentedl200lh . The distribution of 
the sample on the color magnitude diagram (CMD) is plotted in 
Figure [TJ Given the rapid evolution of massive stars off the Main 
Sequence, the position of an A-star on the CMD provides a method 
to estimate the age of the system based on a comparison with theo- 
retical isochrones. The inferred age of the system from the CMD is 
combined with the dynamical system mass from the orbit and sys- 
tem photometry from the literature to test mass-magnitude relations 
at the corresponding age. 
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Figure 1. A colour-magnitude diagram of the 26 stars discussed within this 
work, plotted alongside three different sets of theoretical isochrones at 100, 
500, and 800 Myrs. Those targets with new or refined orbits are plotted 
in blue with a diamond symbol (see Table 0), and those targets for which 
further measurements are required are plotted in red with a cross symbol 
(see Table \5\ and §6.3). Three different sets of theoretical isochrones ar e 
plotted a t ages of 100, 500, an d 800 Myrs; ( top^lLeieune & Schaerei] feOOll) . 
(middle) Marigo et al. (2008), and (bottom) \Siess et alj 12000). 



3 OBSERVATIONS 

With AO systems operating on telescopes ranging in diameter from 
the 3m Shane to the 8m Gemini and VLT, near-infrared images 
were obtained on all targets. In most cases, the filter used for the ob- 
servations was a narrow or broadband filter within the K bandpass, 
though some images were taken within the J and H bandpasses. 
Both the primary and secondary of the pairs were unsaturated in 
the AO images, simplifying the astrometry measurements. Table [2] 
details the key characteristics of the instruments used to acquire the 
new AO observations, with the measured pixel scale and orientation 
for each camera. A subset of the observations were obtained from 
the CFHT and ESO Science Archive Facilities. One measurement 
obtained at the Southern Observatory for Astrophysical Research 
(SOAR) as a part of the VAS T survey, and used wit hin this study, 
has been recently published in Hartko pf et alj J2012h . 
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Table 1. Basic properties of sample members 



HIP 


Name 


HR 


HD 


ADS 


Spectral 


Distance 


V T 


K B 












Type 


(PC) 


(mag) 


(mag) 


Orbit Subsample 


5300 


v Phe 


331 


6767 




A3IV 


57.0±2.0 


5.230±0.001 


4.78±0.02 


9480 


48 Cas 


575 


12111 


1598 


A3IV 


35.3±0.6 


4.533±0.002 


4.08±0.13 a 


11569 


i Cas 


707 


15089 


1860 


A5p 


40.7±1.3 


4.496±0.003 


4.25±0.03 


17954 




1188 


23985 


2799 


A2V+... 


56.5±1.8 


5.259±0.004 


4.81±0.02 


28614 


fj, Ori 


2124 


40932 


4617 




47.5±1.5 


4.150±0.002 


3.64±0.26 


36850 


Castor 


2891/2890 


60178J 


6175 


A2Vm 


15.6±0.9 


1.590±0.020 


1.47±0.03 a 


44127 


lUMa 


3569 


76644 




A7V 


14.51±0.03 


3.159±0.002 


2.67±0.03 b 


47479 




3863 


84121 




A3IV 


72.6±2.2 


5.334±0.003 


4.80±0.02 


76952 


7&B 


5849 


140436 


9757 


B9IV+... 


44.8±1.0 


3.819±0.002 


3.67±0.23 


77660 


bSer 


5895 


141851 




A3Vn 


49.8±0.8 


5.112±0.002 


4.70±0.02 


80628 


v Oph 


6129 


148367 




A3m 


41.0±1.5 


4.657±0.003 


4.17±0.04 


82321 


52 Her 


6254 


152107 


10227 


A2Vspe... 


55.3±1.0 


4.833±0.003 


4.57±0.02 


93506 


CSgr 


7194 


176687 


11950 


A2.5Va 


27.0±0.6 


^\ /Tin i r\ r\f\~) 

2.617±0.003 


2.29±0.23 










Monitorin; 


> Subsample 








128 






224890 




Am... 


70.8±1.7 


6.521±0.004 


6.02±0.02 


2381 




118 


2696 




A3V 


53.1±0.8 


5.188±0.003 


4.83±0.02 


2852 


BGCet 


151 


3326 




A5m... 


48.9±0.8 


6.093±0.004 


5.42±0.01 


5310 


79 Psc 


328 


6695 




A3V 


47.3±1.7 


5.581±0.003 


5.22±0.02 


18217 




1192 


24141 




A5m 


50.5±1.1 


5.806±0.004 


5.37±0.02 


29852 




2265 


43940 




A2V 


61.9±1.0 


5.895±0.003 


5.44±0.02 


51384 




4062 


89571 




F0IV 


40.6±0.6 


5.549±0.003 


4.85±0.02 


65241 


64Vir 


5040 


116235 




A2m 


65.9±1.4 


5.897±0.003 


5.62±0.02 


66223 




5108 


118156 


8956 


F0IV 


69.8±2.0 


6.394±0.004 


5.88±0.02 


103298 


16 Del 


8012 


199254 


14429 


A4V 


60.5±1.1 


5.555±0.003 


5.19±0.02 


109667 




8464 


210739 




A3V 


63.5±2.1 


6.209±0.005 


5.74±0.02 


110787 


p 1 Cep 


8578 


213403 




A2m 


63.2±0.9 


5.857±0.004 


5.54±0.03 


116611 


75 Peg 


8963 


222133 




AlVn 


71.4±1.4 


5.483±0.003 


5.42±0.02 



a - K"-band photometry from Ducati (2002) 

b - X-band photometry from Morel & Magnenat (1978) 




0.1 1.0 10.0 

Separation (arc sec) 



Figure 2. The magnitude difference between primary and secondary for 
each binary system within this study as a function of angular separation. 
The targets are divided into two distinct subsamples, those with new or re- 
fined orbits estimated within this study (filled circles), and those for which 
further measurements are required before an orbital determination can be 
attempted (filled stars). The majority of the systems within the second sub- 
sample were newly resolved as a part of the VAST survey, demonstrating 
the higher dynamic range possible with AO imaging. The filled circles con- 
nected with a solid line represent companions within the same hierarchical 
triple (HIP 44127 and HIP 82321). 



4 DATA ANALYSIS 
4.1 AO image processing 

The AO science images obtained were processed with standard im- 
age reduction steps including dark subtraction, flat fielding, inter- 
polation over bad pixels and sky subtraction. To align all the im- 
ages, the centroid of the bright primary was obtained in each ex- 
posure by fitting a Gaussian to the core of the central point spread 
function. For each system resolved within the observations, an em- 
pirical PSF was determined from the radial profile of the primary, 
after masking any close companion. The empirical PSF was then 
fit to the position and intensity of both components of the system, 
providing a measure of the separation, position angle, and magni- 
tude difference. Uncertainties within the photometry and astrome- 
try were estimated from the standard deviations of the photometric 
and astrometric measurements from each individual exposure be- 
fore combination. 

To ensure accurate determination of the separation and po- 
sition angle, the pixel scale and orientation of the detectors were 
calibrated based on observations of the Trapezium cluster, with 
the exception of data obtained with IRCAL and PHARO. De- 
pending on the total field-of-view of the detector, the positions 
of 20 to 40 Trape zium members were compared w ith the coor- 
dinates reported in iMcCaughrean & Staufferl <fl994h . The average 
derived pixel scale and orientation were computed, and the stan- 
dard deviation of these values were used as the associated errors; 
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Table 2. Observing run details with calculated astrometric calibration values 



Telescope 
(Instrument) 


Date 


Filter 




Plate Scale 
(mas/px) 


True North 

(°) 


Proposal ID 


CFHT 


05/05/2001 -07/05/2001 


H continuum 


34.87±0.14 


1.47±0.19 


2001AF11 


(KIR) 


06/05/2001 


J continuum 


34.82±0.15 


1.51±0.17 


2001 AF11 




13/06/2008 


H 2 (t) = 1 - 


0) 


34.81±0.10 


-2.42±0.13 


2008AC22 




14/06/2008 


Fell 




34.76±0.08 


-2.40±0.11 


2008AC22 




30/08/2009 - 01/09/2009 


H 2 (v = 1 - 


0) 


34.78±0.10 


-2.40±0.20 


2009BC06 




04/02/2010-05/02/2010 


H 2 (v = 1 - 


0) 


34.75±0.08 


-2.37±0.11 


2010AC14 




2011 


H 2 (« = 1 - 


0) 


34.77±0.06 


-2.35±0.09 


2011AC11 


Gemini N 


08/09/2008- 14/11/2008 


Br7 




21.28±0.12 


0.52±0.29 


GN-2008A-Q-74 


(NIRI) 


19/12/2009 


Br7 




21.28±0.12 


0.52±0.29 


GN-2008B-Q-119 




25/06/2010-26/08/2010 


Br7 




21.28±0.12 


0.52±0.29 


GN-2010A-Q-75 


Lick 


24/07/2008 - 25/07/2009 


Br7 




74.55±1.00 


1.68±0.50 


_ 


(IRCAL) 


16/10/2008 - 17/10/2008 


Br7 




74.55±1.00 


1.68±0.50 


„ 


VLT-UT4 


30/06/2004 


IB218 




27.02±0.14 


0.24±0.30 


272.D-5068(A) 


(NACO) 


10/01/2005 


IB218 




27.02±0.15 


0.26±0.30 


074.D-0 180(A) 




08/02/2005 


K s 




27.02±0.15 


0.26±0.30 


074.D-0 180(A) 




08/11/2005 


K s 




26.98±0.08 


0.21±0.20 


076.C-0270(A) 




10/1 1/2005 


K s 




26.98±0.08 


0.21±0.20 


076.D-0108(A) 




06/12/2005 


IB218 




26.98±0.08 


0.21±0.20 


076.D-0108(A) 




04/01/2006 


K s 




27.01±0.13 


0.31 ±0.29 


076.D-0 108(A) 




07/01/2006 


IB218 




27.01±0.13 


0.31 ±0.29 


076.D-0 108(A) 




27/04/2006 


IB218 




27.00±0.12 


0.09±0.26 


077.D-0147(A) 




20/09/2007 


IB218 




26.99±0.12 


-0.23±0.21 


080.D-0348(A) 




20/09/2007 


Ks 




26.99±0.12 


-0.23±0.21 


080.D-0348(A) 




26/09/2007 


IB218 




26.99±0.12 


-0.23±0.21 


080.D-0348(A) 




04/11/2007 


K s 




27.00±0.10 


-0.18±0.20 


080.D-0348(A) 




15/02/2008 


IB218 




27.01±0.14 


-0.12±0.29 


080.D-0348(A) 




24/02/2008 


K s 




27.01±0.14 


-0.12±0.29 


080.D-0348(A) 


Palomar 


1 1/04/2008 


CH 4 




25.00±1.00 


0.70±1.00 




(PHARO) 


12/04/2008 


Br7 




25.00±1.00 


0.70±1.00 






12/07/2008 


CH 4 




25.00±1.00 


0.70±1.00 






13/07/2008 


Br7 




25.00±1.00 


0.70±1.00 





KIR - (Dovon et al. 1998) 

NIRI (Near InfraRed Imager and Spectrometer) - (Hodano et al. 2003) 
IRCAL - (Lloyd et al. 2000) 

NACO (Nasmyth Adaptive Optics System/ Near-Infrared Imager and Spectrograph) - (Lenzen et al. 2003; Rousset et al. 2003) 
PHARO (Palomar High Angular Resolution Observer) - (Hayward e£aL2001) 



the results are given in Table [2] For the data obtained with IR- 
CAL and PHARO, the pixel scale and orientation were calibrated 
from binary systems also observed with instruments calibrated with 
Trapezium measurements. 

4.2 Orbital determination 

For the 13 binaries with sufficient coverage of the orbit, a fit was 
performed for the orbital elements and an estimate of the dynam- 
ical mass was determined. The measurements presented within 
this study were combined with previous measurements con tained 
within the Washington Double Star (WDS: lMason et al.l20"oih Cat- 
alog. These archive measurements were obtained using a variety 
of observational techniques, and date back to the 18 th Century. 
As in some cases the statistical uncertainties were not provided in 
the WDS Catalog, we searched the literature for the formal errors 
for each individual measurement, and only separation and posi- 
tion angle values for which uncertainties could be assigned were 
included within the fitting procedure. A detailed listing of the in- 
dividual measurements used for the orbital determination will be 
made available at the Strasbourg astronomical Data Center (CDS - 
http://cds. u-strasbg . f r). 



Our orbit fitting approach utilises the method presented by 
Hilditcfl l l200lh . and demonstrated b y an applica t ion to mea- 
surements of the T Tau S system jKbhler et alj [2008). This 
method is simi l ar to t he grid-based search technique developed by 
Hartk opf et al.1 (Q~989). At each epoch of observation ti, the Xi,yi 
position of the companion with respect to the primary is measured 
in the observed tangent plane. These values are related to the true 
position of the secondary in the orbital plane (x'^y'i) through the 
following equations 

Xi = Ax[ + Fy'i 
Vi = Bx[ + Gy'i 

where A, B, F and G are the orbital Thiele-Innes elements, with 

A = a(cos uj cos f2 — sin uj sin il cos i) 
B = a (cos a; sin + sin a; cos cos i) 

(2) 

F = a(— sin uj cos Q. — cos uj sin Q cos i) 
G = a(— sin a; sin Q + cos uj cos Q cos i) 

where a is the semi-major axis of the orbit, uj the longitude of peri- 
astron, Q, the longitude of the ascending node, and i the inclination 
- four of the seven orbital elements. The position of the compan- 
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ion in the orbital plane (x^y'i) can also be expressed through the 
remaining orbital elements (e, P, To) as 



Table 3. Summary of theoretical model grids 



x'j — cos E — e 

y[ — \/l — e 2 sin E 



(3) 



where e is the eccentricity of the system. The eccentric anomaly 
(E) can be determined from a numerical solution to Kepler's equa- 
tion 



M = E - e sin E 

= (2n/P)(U-T Q ) 



(4) 



where M is the mean anomaly, P the period of the system, and To 
the epoch of periastron passage. At each epoch of observation, the 
position of the component in the orbital plane can be defined using 
just three of the orbital elements (e, P, and To). The orbital posi- 
tion at each epoch can then be converted into the observed position 
using the equations in Equation [TJ through a least-squares determi- 
nation of the four Thiele-Innes elements. 

An initial estimation of the orbital parameters of each sys- 
tem can be determined through an iterative three-dimensional grid 
search of e, P, and To. A wider range of parameter values were 
searched, with 100 linear steps searched over a range of ^ e < 1, 
500 linear steps between 1 ^ \og(P[yrs]) ^ 3, and To initially 
distributed between the years 2000.0- (P/2) and 2000.0+ (P/2). 
At each position within this three-dimensional grid, the fit orbital 
positions (x'^y'i) were directly calculated (Eqn. O, with the four 
remaining orbital parameters (a, i, uj, Q) estimated from a least- 
squares fit to the observed positions using Equations[T]and|2] After 
the \ 2 statistic was calculated at each position within the grid, the 
range of To values searched was reduced by a factor of 10 centred 
on the optimum value of To found within the previous iteration. 
This process was repeated until the step size in To was reduced to 
less than one day. The values fo r a, i, uj and fl can be determined 
from an inversion of Equation |2]( Green 1985]), 



arctan 



B 



uj = ^arctan 



2 arctan 



G + A 
B -F 
G + A 



arctan 



B + F 

G-A 



(B + F) sin (w + O) 



B 



(B 
F 



F) sin (u) - Q) 



(5) 



sin (uj + fi) (1 + cos (i)) 

The orbital parameters calculated at each position within the (P, e) 
grid are then used as a starting point for a Levenberg-Marquardt 
minimisation to ensure the minimum of the \ 2 distribution is 
found. The set of orbital parameters with the minimum \ 2 statistic 
was then used as the orbit solution for the system. 

The shape of the x distribution in the vicinity of the global 
minimu m can be used to d etermine the la uncertainties of each pa- 
rameter jPress et al.|[l992h . By perturbing an individual parameter 
away from the global minimum, and optimising the remaining pa- 
rameters, a region of the x 2 distribution can be calculated where the 
X 2 statistic is less than Xmin + !• This region encloses 68% of the 
probability distribution, and is not necessarily symmetric about the 
minimum x 2 value. Our implementation of t he orbit fitting method 
was tested against four well studied systems (Bonnefo v et alj2009l : 
IDuduv et alJl2009l : iLiu et al.ll2008h . with the resulting parameters 
being within the published la uncertainties. 



Grid 
reference 



Mass Range 
(M ) 



Metallicity 

(Z) 



Lejeune & Schaerer (2001) 


0.80 


< 


M 


< 


5.0 


0.02 


Mariao et al. (2008) 


0.15 


< 


M 


< 


5.0 


0.02 


Siess et al. (2000) 


0.10 


< 


M 


< 


5.0 


0.02 


Baraffe etal. (1998) 


0.08 


< 


M 


< 


1.5 


0.02 



4.3 Theoretical mass-magnitude relations 

We obtained four different grids of evolutionary models with which 
the dynamical system masses estimated from the fitted orbital 
parameters were compared. Table [3] lists the mass range, metal- 
licity, and literature reference for each of the four grids. The 
iLeieune & Schaereri feOOll) . iMarigo et af] (2008), and ISiess et al.1 
(2000) grids covers a significant portion of the lifespan of a typ- 
ical A-type star, and as such we applied a ma ximum age cut-off a t 
1 Gyr. In addition to these grids, models from lBaraffe et alj i 19981) 
were obtained in order to study a pair of lower-mass companions 
presented in §6.1.2. Each grid was converted into the photomet- 
ric systems used within this study - Tycho V and 2MASS Ks 
(Carpenter 2001), before producing a high resolution (dM/M© = 
0.001) mass-magnitude relation, created through cubic interpola- 
tion of the grid data, as shown in Figure [5] 

The absolute V- and if-band magnitudes were calculated 
from the V-band magnitude differences obtained from the litera- 
ture, and the if -band magnitude differences presented within this 
study (Tables [4] and |5j. The individual component magnitudes and 
V — K colour for each system are presented in Table [6] Within 
the A-type star mass range, the mass-magnitude relations signif- 
icantly change as a function of the age of the system due to the 
rapid evolution of A-type stars across the CMD. The age of each 
system is therefore estimated, based on the position of the primary 
on the CMD (Figure [4}, before a mass for each component is es- 
timated from the mass-magnitude relations. The estimated masses 
of each component are summed to produce an estimation of the 
system mass, hereafter called the theoretical system mass. 

To demonstrate the analysis procedure, a hypothetical binary 
system of magnitudes Mk(A) = 1.25, Mk(B) = 5.50 was 
used to construct the mass-age relation for each component, and 
their corresponding sum (Figure |5j- The evolution of the mass- 
magnitude relation as a function of age can be then visualised as 
a continuous function for both components within the system. For 
this example a set of models were used which includes the pre- 
Main Sequence (PMS) contraction phase of lower-mass stars, as 
demonstrated by the increase in the derived mass as a function 
of age for the companion. The two mass-age curves can then be 
summed to produce a mass-age relation for the system in question. 
Using the age estimated for the system based upon its position on 
the CMD, a theoretical system mass can be estimated and com- 
pared with the dynamical system mass determined from the orbital 
elements. 



5 RESULTS 

5.1 Astrometric results 

The astrometry and photometry measurements of the two subsam- 
ples are given in Tables|4]and[5] Both tables contain the Hipparcos 
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Table 4. Measured binary position angle and separation for previously resolved systems 



HIP 


Comp. 


WDS 


Discoverer 


Previous 


Instrument 


Epoch 


e 


p 


AK 






Designation 


Designation 


Orbit 






(deg) 


(arcsec) 


(mag) 


5300 


AB 


01078-4129 


RST3352 


Soderhielm (1999) 


NACO 


2005.02 


15.73±0.76 


0.089±0.001 


0.5()7±0.020 












NACO 


2007.72 


126.69±0.38 


0.114±0.001 


*0.573±0.018 


9480 


A-B 


02020+7054 


BU513 


Mason et al. (19991 


IRCAL 


2008.79 


293.97±0.75 


0.679±0.011 


1.110±0.056 










KIR 


2009.67 


297.36±0.22 


0.666±0.002 


1.172±0.013 












NIRI 


2010.65 


302. 11 ±0.29 


0.644±0.004 


*1.137±0.007 


11569 


Aa,Ab 


02291+6724 


CHR6 


Drummond et al. (2003) 


IRCAL 


2008.79 


43.50±0.80 


0.575±0.006 


*1.906±0.026 










KIR 


2010.10 


41.52±2.61 


0.602±0.023 


1.771±0.301 


17954 


AB 


03503+2535 


STT65 


Docobo & Ling (2007) 


NIRI 


2008.87 


194.76±0.31 


0.200±0.001 


*0.232±0.005 


28614 


AB 


06024+0939 


A 2715 


Muterspaugh et al. (2008) 


NIRI 


2009.96 


21.95±0.29 


0.403±0.002 


*1.008±0.006 










KIR 


2011.31 


21.31±1.06 


0.397±0.006 


0.896±0.120 


36850 


AB 


07346+3153 


STF1110 


Docobo & Costa (1985) 


KTR 


2010.09 


57.46±0.14 


4.663±0.013 


*0.698±0.050 


44127 


BC 


08592+4803 


HU628 


Eggen (1967) 


PHARO 


2008.28 


242.01 ±1.04 


0.433±0.017 


*0.005±0.031 










KIR 


2010.10 


228.98±0.27 


0.568±0.003 


0.027±0.023 












KIR 


2011.30 


223.37±0.31 


0.656±0.003 


-0.011±0.039 


47479 


AB 


09407-5759 


B 780 


Soderhielm (1999) 


NACO 


2008.12 


300.22±0.58 


0.099±0.001 


*0.082±0.017 


76952 


AB 


15427+2618 


STF1967 


Muterspaugh et al. (2010) 


PHARO 


2008.28 


112.65±1.07 


0.683±0.028 


_ 












KIR 


2011.31 


111.37±0.33 


0.620±0.005 


*1.056±0.030 


77660 


AB 


15513-0305 


CHR51 


Docobo etal. (2010) 


KIR 


2001.34 


82.72±2.42 


0.188±0.018 


_ 










KIR 


2001.34 


81.97±1.59 


0.197±0.008 


- 












KIR 


2001.35 


83.38±3.67 


0.185±0.O16 














NACO 


2004.50 


71.87±0.31 


0.247±0.001 


*1.487±0.006 












NACO 


2006.32 


69.74±0.42 


0.274±0.003 


1.472±0.097 












KIR 


2010.10 


66.69±0.47 


0.330±0.004 


1.483±0.058 












KIR 


2011.30 


64.21±2.57 


0.340±0.014 


1.455±0.194 


80628 


Ab,Ab 


16278-0822 


RST3949 




PHARO 


2008.28 


22.56±1.05 


0.674±0.027 


*2.072±0.066 












IRCAL 


2008.56 


23.53±0.71 


0.685±0.009 


1.978±0.057 












KIR 


2011.31 


34.37±0.29 


0.771 ±0.003 


2.237±0.038 


82321 


BC 


16492+4559 


A 1866 


Popovic (1969) 


PHARO 


2008.53 


247.33±1.03 


0.286±0.012 












KIR 


2011.30 


260.51 ±0.53 


0.301 ±0.005 


*0.299±0.044 


93506 


AB 


19026-2953 


HDO 150 


Mason et al. (1999) 


KIR 


2008.45 


31.05±2.31 


0.194±0.010 














KIR 


2011.31 


285.59±0.71 


0.304±0.004 


*0.320±0.043 



* - Photometry measurement used to determine the magnitude of each component in Table[6] 



Table 5. Measured binary position angle and separation for potential orbital monitoring targets 



HIP 


Instrument 


Epoch 


e 


P 


AK 


Projected 








(deg) 


(arcsec) 


(mag) 


Separation (AU) 


128 


NIRI 


2008.72 


80.64±0.29 


0.983±0.006 


*3.521±0.021 


69.59±1.72 




IRCAL 


2008.79 


80.73±0.52 


0.975±0.013 


3.355±0.037 


68.98±1.90 


2381 


NACO 


2005.93 


279.51 ±0.26 


1.765±0.006 


*6.206±0.088 


93.75±1.48 




NACO 


2007.73 


279.13±0.22 


1.767±0.008 


6.229±0.055 


93.86±1.51 




NIRI 


2008.79 


278.54±0.29 


1.759±0.()13 


7.401±0.711 


93.39±1.59 


2852 


NIRI 


2008.79 


260.55±0.30 


0.931±0.006 


*5.001±0.056 


45.51±0.82 


5310 


NIRI 


2008.79 


174.99±0.30 


0.357±0.002 


*3.7I9±0.023 


16.88±0.62 


18217 


IRCAL 


2008.79 


64.78±0.65 


1.037±0.012 


2.301±0.029 


52.37±1.24 




NIRI 


2008.87 


64.96±0.29 


1.033±0.006 


*2.41±0.010 


52.13±1.12 


29852 


NACO 


2005.85 


210.76±0.22 


0.218±0.001 


*2.003±0.013 


13.49±0.21 




NACO 


2005.86 


210.80±0.21 


0.217±0.001 


1.983±0.012 


13.44±0.21 


51384 


PHARO 


2008.18 


212.43±1.05 


2.078±0.084 


*4.543±0.169 


84.44±3.66 


65241 


NACO 


2005.10 


196.95±0.62 


0.328±0.003 


3.061±0.231 


21.65±0.50 




NACO 


2008.15 


212.75±0.57 


0.258±0.003 


*3.025±0.097 


17.01±0.42 


66223 


PHARO 


2008.53 


187.71±1.02 


1.381±0.056 


*5.660±0.105 


96.42±4.73 


103298 


NIRI 


2008.69 


115.72±0.30 


0.220±0.001 


*2.962±0.007 


13.29±0.25 


109667 


NIRI 


2008.69 


285.17±0.30 


1.117±0.007 


*3.997±0.011 


70.88±2.38 




KIR 


2009.66 


284.66±0.25 


1.108±0.006 


4.051 ±0.073 


70.28±2.35 




NIRI 


2010.48 


284.38±0.29 


1.101±0.007 


4.044±0.103 


69.84±2.34 


110787 


NIRI 


2008.71 


211.11±0.35 


0.291±0.002 


*3.891±0.043 


18.38±0.29 


116611 


NIRI 


2008.75 


173.11±0.32 


0.950±0.006 


*5.930±0.089 


67.78±1.38 




NIRI 


2010.48 


172.11±0.34 


0.943±0.006 


5.794±0.076 


67.31±1.34 



* - Photometry measurement used to determine the magnitude of each component in Table[6] 
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Table 6. Derived parameters of individual components 



HIP 


Comp. 


V T 




V — K 


M v 


M K 


Additional 


SB 


Reference 






(mag) 


(mag) 




(mag) 


(mag) 


Components 


Type 




Orbit Subsample 


5300 


A 


5.51±0.02 


5.28±0.02 


0.23±0.03 


1.73±0.08 


1.50±0.08 










B 


6.85±0.08 


5.85±0.02 


0.99±0.08 


3.07±0.11 


2.08±0.08 








9480 1 " 


A 


4.69±0.01 


4.40±0.13 


0.29±0.13 


1.95±0.03 


1.67±0.13 










B 


6.72±0.01 


5.54±0.13 


1.18±0.13 


3.98±0.04 


2.80±0.13 








11569 t 


Aa 


4.67±0.01 


4.79±0.03 


-0.12±0.03 


1.62±0.07 


1.74±0.08 










Ab 


8.65±0.16 


6.71 ±0.03 


1.94±0.16 


5.60±0.17 


3.66±0.08 








17954 


A 


5.81±0.04 


5.46±0.02 


0.35±0.04 


2.05±0.08 


1.70±0.07 










B 


6.26±0.06 


5.69±0.02 


0.57±0.06 


2.50±0.09 


1.93±0.07 








28614 


A 


4.32±0.01 


4.00±0.26 


0.33±0.26 


0.94±0.07 


0.61±0.27 


Aa, Ab 


SB1 


Fekel et al. (2002) 




B 


6.22±0.02 


5.01±0.26 


1.22±0.26 


2.84±0.07 


1.62±0.27 


Ba, Bb 


SB2 


Fekel et al. (2002) 


36850 


A 


1.98±0.02 


1.93±0.03 


0.05±0.04 


1.02±0.13 


0.97±0.13 


Aa, Ab 


SB1 


Vinter Hansen et al. (1940) 




B 


2.88±0.03 


2.63±0.04 


0.25±0.05 


1.92±0.13 


1.66±0.13 


Ba, Bb 


SB1 


Vinter Hansen et al. (1940) 


44127 


A 


3.16±0.02 


2.71 ±0.03 


0.45±0.04 


2.35±0.02 


1.90±0.03 








B 


10.88±0.11 


6.90±0.04 


3.98±0.12 


10.07±0.11 


6.09±0.04 










C 


11.08±0.12 


6.90±0.04 


4.18±0.12 


10.27±0.12 


6.09±0.04 








47479 


A 


5.82±0.01 


5.51±0.02 


0.31 ±0.02 


1.51 ±0.07 


1.21 ±0.07 










B 


6.45±0.01 


5.59±0.02 


0.85±0.03 


2.14±0.07 


1.29±0.07 








76952 


A 


4.05±0.01 


4.02±0.23 


0.03±0.23 


0.80±0.05 


0.76±0.23 










B 


5.60±0.02 


5.07±0.23 


0.53±0.23 


2.35±0.05 


1.82±0.23 








77660 


A 


5.21±0.01 


4.94±0.02 


0.26±0.02 


1.72±0.04 


1.46±0.04 










B 


7.82±0.04 


6.43±0.02 


1.38±0.04 


4.33±0.05 


2.95±0.04 








80628 


Aa 


4.68±0.01 


4.32±0.04 


0.37±0.04 


1.62±0.08 


1.25±0.09 


Aal, Aa2 


SB2 


Gutmann (1965) 




Ab 


8.80±0.10 


6.39±0.07 


2.41±0.12 


5.74±0.13 


3.33±0.11 






82321 


A 


4.87±0.01 


4.73±0.02 


0.14±0.02 


1.16±0.04 


1.02±0.04 










B 


9.19±0.11 


7.33±0.03 


1.86±0.11 


5.47±0.12 


3.62±0.05 










C 


9.29±0.11 


7.63±0.04 


1.66±0.12 


5.57±0.12 


3.92±0.05 








93506 


A 


3.27±0.01 


2.90±0.23 


0.37±0.23 


1.11 ±0.05 


0.74±0.24 










B 


3.48±0.01 


3.22±0.24 


0.26±0.24 


1.32±0.05 


1.06±0.24 








Monitoring Subsample 


128 


AB 




6.06±0.02 






1.81 ±0.06 


A, B 


SB1 


Carauillat et al. (2003) 




C 




9.58±0.03 






5.33±0.06 








2381 


A 




4.83±0.02 






1.21±0.04 










B 




11.04±0.09 






7.41±0.10 








2852 


A 




5.43±0.02 






1.98±0.04 










B 




10.43±0.06 






6.98±0.07 








5310 


A 




5.25±0.02 






1.88±0.08 










B 




8.97±0.03 






5.60±0.08 








18217 


A 




5.48±0.02 






1.97±0.05 










B 




7.90±0.02 






4.38±0.05 








29852 


A 




5.60±0.02 






1.64±0.04 










B 




7.60±0.02 






3.64±0.04 








51384 


A 




4.87±0.02 






1.82±0.04 










B 




9.41±0.17 






6.37±0.17 








65241 


A 




5.69±0.02 






1.59±0.05 










B 




8.71 ±0.09 






4.62±0.11 








66223 


Aa 




5.89±0.02 






1.67±0.06 










Ab 




11.55±0.11 






7.33±0.12 








103298 


Aa 




5.26±0.02 






1.35±0.04 










Ab 




8.22±0.02 






4.31 ±0.04 








109667 


A 




5.76±0.02 






1.75±0.07 










B 




9.76±0.02 






5.75±0.07 








110787 


A 




5.57±0.03 






1.57±0.04 










B 




9.46±0.05 






5.46±0.06 








116611 


Aa 




5.42±0.02 






1.16±0.05 


Aal, Aa2 


SB1 


Rucinski et al. (2005) 




Ab 




11.35±0.09 






7.09±0.10 







t - The spectroscopic binary reported by Abt (1965) is the system resolved within the AO data. 

t - The individual component magnitudes may be significantly biased due to the presence of additional companions within the resolution limit of the 
Tycho2 and 2MASS observations. 
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Figure 3. Mass-magnitude relations were constructed from each of the four 
model grids. The evolution of the mass-magnitude relation is shown for 
each grid, with the 10 Myr (black line), 100 Myr (blue line) and 1 Gyr (red 
line) relations plotted. For reference the dashed line indicates the faintest 
A-type star within the sample, and the dotted line indicates the Zero-Age 
Main Sequence magnitude of an M0 star. A-type stars typically have an 
absolute X-band magnitude ranging between Mk = and Mjf = 2, 
wh ere the mass is sign i ficantl y de pendent on the age of the star. Within 
the lLeieune & SchaereJ 1200 ll) and lMarigo et alj 12008) models, the mass- 
magnitude relation of stars fainter tha n Mj< ~ 2 is not depe ndent on the 
age of the star. The lSiess et al] l2000l) and lBaraffe et alj Il998h models in- 
clude a description of the contraction phase of lower-mass stars onto the 
Main Sequence during the early portion of its life. 

designation of the primary, the components of the system under in- 
vestigation, the instrument and epoch of observation, and the mea- 
sured astrometric values with corresponding uncertainties. For the 
orbit subsample, the WDS designation and discoverer code are also 
listed for reference. 

Combining the A'-band magnitudes of the sample with the 
A A values reported in Tables [4] and [5] and the Hipparcos paral- 
lax, allowed for an estimation of the A'-band apparent and abso- 
lute magnitudes of the resolved components (Table [6). For sys- 



Figure 4. The age of each system can be estimated based on the 
position of each c o mpone nt on the c olour-magnitu d e diag ram (top - 
Leieune & SchaereJ 1200 ll) . middle - iMarigo et all 12003) . bottom - 
Siess et alj (2000)). This procedure, while demonstrated for only two sys- 
tems in this figure, was repeated for all the remaining members of the orbit 
subsample. For each system, the age estimated from each of the three model 
grids is reported in Table llOl 



terns with AV measurements reported within the literature (Ta- 
ble [7J, the corresponding estimated V-band apparent and absolute 
magnitudes, and V — K colours for each resolved component are 
reported. Seven members of the overall sample are hierarchical 
systems with at least one of the components resolved within our 
dataset consisting of multiple sub-components, indicated in Table 
[6] An example of this is HIP 128; our AO data are able to resolve 
a previously-unknown binary companion within this system (HIP 
128 C at ~ l'.'O), but are of insufficient angular resolution to re- 
solve the previously-known spectroscopic component HIP 128 B. 
Without an estimate of the AV' or AK between HIP 128 A and HIP 
128 B, the individual magnitudes cannot be estimated and therefore 
only the blended magnitudes of the two components are listed. 
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Table 8. Estimated orbital parameters and corresponding la uncertainties 



HIP 


N. 


Period 


Semi-major 


Inclination 


Longitude 


Epoch of 


Eccentricity 


Longitude 


System 




Meas. 




axis 




of Node 


Periastron 




of Periastron 


Mass 






f (,yrs; 


n ("1 


« *. ) 


\i { ) 






LU ( ) 


U vi O ) 


5300 


19 


28.36±g°| 


2396+ - 0005 
v.£jyv_Q qqq5 


65 - 3 ±03 


142.718:1 


2007.94+ 8'8| 


474+0.002 
— 002 


333.9l8'| 


3.16±0.34 


9480 


66 




n £.4+0:002 
U.014_ Q02 


16 7+° :9 


48.2±|f 


1964.35+8-09 


isc+oiooi 


19.5li;s 


2.72±0.13 


11569 


11 


5o.2l i;8 


429+0.007 


149.0+ \'l 


180.0l| : | 


1993.2418:81 


642+0- 009 


331.3l|-f 


2.12±0.25 


17954 


44 


61.2l8;i 


447+0.002 
u -™-0 .002 




26 4+° :2 
zo -^-0.2 


1998.3018:81 


n fi7S+0.002 
U-628_ 002 


340.3+8 : ? 


4.15±0.39 


28614 


69 


18 A4 1+0.008 
18.641_ 0Qg 


2741 + ' 0001 
u.z/t i_ 0Q01 


96.59l°;8| 


74 7fi+0.02 


2003.74218 88! 


O744+O.OOI 
u - '^-0.001 


217.08±8;8I 


6.36±0.62 


36850 


207 


466.8 + |? 


a 7 S +0.05 


113.56±8;8i 


41.2l8 : l 


1957.318:3 


13Q + 0.007 
u-J->J_0.006 


249.318:1 


5.42±0.97 


44127 


9 


39.4±8i 


ri 70+O.OI 
u/u -0.01 


111.61H 


24.51H 


1999.1 + 8 ? 


nz+0.02 
u --"-0.02 


354.2+lf 


0.68±0.04 


47479 


5 


10.74±0;8t 


n 1?n7+0.0006 
u.izu/_ 0006 


128.41 


87 4+0.5 


2005.8918^81 


o %5+o-o° 9 


18.6ll:i 


5.83±0.53 


76952 


104 


ni 2+0-4 


n 79Q+0.006 
"•' zy _0.006 


9 44 c+0:07 
"•'•■-J -0.07 


iii-75l8:8I 


1931.618:1 


0.4818:8! 


io3.8l8:f 


4.19±0.30 


80628 


8 


82.81 13 


79+0.03 


3i-2ll :u 




1994.llo : 9 


0-4518:81 


177.9+6-f, 


4.99±0.75 


82321 


12 




779+0.005 


37.41 1:| 


57.5l| : ? 


1991.2l2 : f 


13+0.02 
u - 1J -0.02 


69.41 f° 


1.16±0.09 


93506 


42 




O4S9+O.OOI 


iii.il8:l 
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Figure 5. The mass-magnitude relation for a hypothetical binary system 
with M K (A) = 1.25 ± 0.05 and M K (B) = 5.50 ± 0.10, as a function 
of system age. The evolution of the mass-magnitude relation is shown for 
both components (A - blue region, B - red region), and the system as a 
whole (grey region). The extent of the region in each case represents the 
uncertainty in the mass estimate due to the uncertainties of the magnitudes 
of each component. 



5.2 Orbital elements and dynamical masses 

The orbital parameters for each system with sufficient orbital cov- 
erage are listed in Table[8] alongside a system mass estimated from 
Kepler's Third Law (hereafter the dynamical system mass) and the 
number of measurements used to fit the orbit. In Figure [6] the re- 
fined orbital fits incorporating the new data are plotted along with 
the previously reported orbit (references listed within Tabled- The 
resulting orbits span a range of periods from 10.78 to 467.4 years, 
and a range of semi-major axes between (/.'12 to 6'.'78. The rela- 
tively short period of each system allows for orbital motion to be 
resolved over very short baselines, typically on the order of months. 
The changes to the estimated dynamical system mass between the 
previously published orbit fit and the refined fit presented within 
this study are shown in Table [9] 



Table 7. Visual magnitude differences for a subset of the sample 



HIP 


AV 


Reference 


5300 


1.34±0.10* 


Horchetal. (2001) 


9480 


2.03±0.01 


Fabricius & Makarov (2000) 


11569 


3.98±0.16 


Christou & Drummond (2006) 


17954 


0.45±0.10* 


Horch et al. (2004) 


28614 


1.90±0.02 


Fabricius & Makarov (2000) 


36850 


0.90±0.03 


Worlev (1969) 


44127 A,BC 


7.06±0.10* 


Baize & Petit (1989) 


44127 BC 


0.20±0.10* 


Mason et al. (2001) 


47479 


0.63±0.02 


Mason et al. (2001) 


76952 


1.55±0.02 


Fabricius & Makarov (2000) 


77660 


2.61±0.04 


Docoboetal. (2010) 


80628 


4.12±0.10* 


Mason et al. (2001) 


82321 A,BC 


3.61±0.10* 


Mason et al. (2001) 


83231 BC 


0.10±0.10* 


Mason et al. (2001) 


93506 


0.21±0.02 


Fabricius & Makarov (2000) 



* - For AV measurements without uncertainties, 0.10 is used 



6 DISCUSSION 

We have presented high resolution observations obtained for 26 
nearby multiple systems with A-type primaries with projected sep- 
arations within 100 AU. The subset of 12 targets with orbit fits have 
been further divided into four distinct categories primarily based 
on a comparison between the dynamical mass and the theoretical 
mass estimated from the mass-magnitude relations, as shown in 
Figure |7] T he theore t ical m a ss estimates for each syste m using the 
iLeieune & Schaererl l l200ll) , iMarigo et al] d2008h . and Isiess et al.1 
fcoooh grids are shown in Table [10] Those systems with only 
two known components are discussed in §6.1.1, where the impor- 
tance of metallicity is described. Two hierarchical systems resolved 
within our data are discussed in §6.1.2 and §6.1.3, which allow for 
a comparison to the models within the K- to M-type spectral range. 
For systems with a dynamical mass excess, significantly higher 
than the mass predicted from the theoretical mass-magnitude re- 
lations, the subset with known spectroscopic components is dis- 
cussed in §6.2.1, and we present three systems with evidence sug- 
gestive of an additional unresolved component in §6.2.2. The re- 
maining targets are discussed in the context of continued monitor- 
ing of the orbital motion in §6.3, of these 1 1 are newly resolved as 



10 De Rosa et al. 




■ ■ ■ ■ ■ i 

0.1 0.0 -0.1 -0.2 -0.5 0.6 0.4 02 0.0 -0.2 -04 -0.6 -0.8 

to (oreseeonds) to (orcstconds) 




0.6 0.4 0.2 0.0 -0.2 -04 0.4 0.2 0.0 -0.2 -0.4 -0.6 

to (orcstconds) to (orcstconds) 




-"'r ■ ■ ■ ■ ■ H -«l ... ■ ... ■ ... ■ ... ■ ... ■ ... I 

0.5 0.2 0.1 0.0 -0.1 -0.2 6 4 2 -2 -4 -6 

to (orcstconds) to (orcstconds) 



Figure 6. Combining our high resolution observations with historical measurements, refined orbits for 6 binary systems are plotted. The previous orbital fit, 
obtained from the Sixth Orbit Catalog, is plotted for reference with a dashed line. Each plot uses a similar symbol scheme to the Sixth Orbit Catalog; open 
blue circles representing eyepiece interferometry, filled blue circles speckle interferometry, green crosses micrometrical observations, and violet asterisks pho- 
tographic measurements. Our high resolution observations presented within this study are plotted as filled red stars. For each measurement, the corresponding 
O-C line is plotted, showing the difference between expected and actual position within the orbital path. Symbols in grey represent those measurements pre- 
sented without formal errors, and are not used while estimating the orbital parameters. Within each plot, the 57mas radius black disc represents the resolution 
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Figure 6. Combining our high resolution observations with historical measurements, refined orbits for 6 binary systems are plotted. The previous orbital fit, 
obtained from the Sixth Orbit Cataloge, is plotted for reference with a dashed line. Each plot uses a similar symbol scheme to the Sixth Orbit Catalog; open 
blue circles representing eyepiece interferometry, filled blue circles speckle interferometry, green crosses micrometrical observations, and violet asterisks pho- 
tographic measurements. Our high resolution observations presented within this study are plotted as filled red stars. For each measurement, the corresponding 
O-C line is plotted, showing the difference between expected and actual position within the orbital path. Symbols in grey represent those measurements pre- 
sented without formal errors, and are not used while estimating the orbital parameters. Within each plot, the 57mas radius black disc represents the resolution 
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Table 10. Comparison of dynamical and theoretical system masses for each system within the orbit subsample 



HIP 


Dynamical 

System 
Mass (Mq) 


Lejeune & Schaerer i M^rrigoetaL 
Age System ■ Age System 
(Myr) Mass(M0) • (Myr) Mass (Mq) 


Siess et al. 
Age System 
(Myr) Mass (Mq) 


Systems identified as doubles (§6.1.1) 


5300 
9480 

11569 AaAbt 
76952 


3.16±0.34 
2.72±0.13 
2.12±0.25 
4.19±0.30 


280 3.56±0.13 400 3.47±0.13 
250 3.15±0.20 280 3.12±0.20 
100 2.45±0.07 100 2.89±0.15 
315 4.38±0.45 400 4.10±0.40 


355 
160 
100 
315 


3.69±0.16 
3.25±0.31 
2.97±0.25 
4.38±0.45 




Low 


mass binary pair within hierarchical system (§6.1.2, §6.1.3) 






44127 BC 
82321 BC 


0.68±0.04 
1.16±0.09 


250 0.88±0.01 
400 1.91±0.02 ■ 445 1.87±0.02 


50 
500 


0.83±0.01 
1.90±0.02 




Syst 


ems with one unresolved spectroscopic component (§6.2.1) 






80628 


4.99±0.75 


630 3.03±0.08 ■ 630 3.06±0.18 


445 


3.15±0.13 


Systems with two unresolved spectroscopic components (§6.2.1) 


28614 
36850 


6.36±0.62 
5.42±0.97 


630 3.98±0.32 ■ 630 3.96±0.28 
280 4.26±0.28 355 4.15±0.24 


500 
315 


4.18±0.34 
4.34±0.31 




Systems 


with evidence of unresolved spectroscopic components (§6.2.2) 






17954 
47479 
93506 


4.15±0.39 
5.83±0.53 
5.26±0.37 


400 3.49±0.11 400 3.44±0.12 
560 3.99±0.10 560 3.96±0.10 
710 4.11±0.25 710 4.07±0.23 


225 
445 
500 


3.67±0.16 
4.17±0.14 
4.45±0.36 



t - The ages and theoretical masses may be significantly biased due to the presence of additional companions 
within the resolution limit of the Tycho2 and 2MASS observations. 



Table 9. Previous and revised dynamical system masses 



28614 h 
4,7479 

t 36850 
93506 - 



17954 



76952 



5300 
'<+< 9480 



i 82321 



1 2 3 4 5 6 
Theoretical Mass (M/M ) 



HIP 


Old Estimate 

(Mq) 


New Estimate 

(Mq) 


5300 


3.45 


3.16±0.34 


9480 


2.98±0.39 


2.72±0.13 


11569 


2.12±0.33 


2.12±0.25 


17954 


3.37±0.35 


4.15±0.39 


28614 


6.33±0.62 


6.36±0.62 


36850 


5.51 


5.42±0.97 


44127 


0.61±0.03 


0.68±0.04 


47479 


7.50 


5.83±0.53 


76952 


4.14±0.28 


4.19±0.30 


80628 




4.99±0.75 


82321 


1.00 


1.16±0.09 


93506 


5.20±0.36 


5.26±0.37 



Figure 7. Using the method introduced in §4.4, a comparison can be made 
between the dynamical mass determined from the orbit and the mass es- 
timated from theoretical mass-mag nitude relations for the stars within the 
orbit subsample (the iMarigo et alj (2008) models are used for this exam- 
ple). The systems with an A-type star primary which are known to consist 
only of two components are denoted as black points (§6.1.1). Two hierar- 
chical systems were fully resolved with our data, and the lower-mass pair of 
each system are in green (§6.1.2). The systems which have a significantly 
discrepant dynamical mass can be explained by the presence of an unre- 
solved spectroscopic companion within our data. The targets with known 
spectroscopic components are plotted in red (§6.2.1), while those systems 
with evidence suggesting a previously unknown spectroscopic component 
are plotted in blue (§6.2.2). The dashed line denotes the equivalence be- 
tween the dynamical and theoretical mass. 



a part of the VAST survey, 2 were resolved in recent multiplic- 
ity surveys of nearby Southern A-type stars jlvanov et alj|2006l : 
lEhrenreich eT^]|2010h. and 1 resolved wi thin a large speckle in- 
terferometry survey " iMcAlisteret all 19871) . 

6.1 Comparison to theoretical models 

6.1.1 A-type binaries 

Four targets within the orbit subsample are systems where the two 
known components have been resolved within our high resolution 
data; HIPs 5300, 9480, 11569, and 76952. For each of these four 
targets, the mass-magnitude relations were used to determine how 
the system mass changes as a function of our estimate of the system 
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Figure 8. The system mass as a function of age based on the mass-magnitude relations derived from three of the model grids. Each panel is similar in nature 
to Figure|5] with the spread in the mass-age relation introduced by uncertainties in both the measured _R"-band magnitude, and the distance determination. The 
horizontal dot-dashed line indicates the dynamical mass determined from the system orbit, with the uncertainty denoted by the line-shaded region. The vertical 
dotted line indicates the age of the system derived from the isochrones, with the line-shaded region denoting the range of ages consistent with the uncertainties 
in the position of the primary on the CMD. The primary of the HIP 1 1569 system is significantly bluer than expected for a Main Sequence star, and as such its 
age has been assigned to 100 Myr, with no corresponding uncertainty. The presence of additional components to the HIP 11569 system, within the resolution 
limit of both the Tycho2 and 2MASS observations, is the likely cause of this bias. 
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Figure 9. (top panel): The dynamical mass of the HIP 5300 system (hor- 
izontal line-shaded region) is significantly below that of the system mass 
estimated from mass-magnitude relati ons (solid shaded region) c alculated 
from the Solar metallicity models from lLeieune & Schaerea i200lh . (middle 
panel): Assuming a higher metallicity (Z = 0.04), the position of the com- 
ponents on the colour-magnitude diagram suggests a younger age (denoted 
by the vertical line-shaded region), and the mass-magnitude relation derived 
from the metal-enhanced models suggests a systematically lower system 
mass, although the dynamical mass is still discrepant, (bottom panel): With 
a metallicity of Z = 0.10, the age estimate is younger still (< 10 Myr), and 
the discrepancy between the dynamical mass and the mass estimated from 
the mass-magnitude relation is removed. 

age, shown in Figure[8] For one system, HIP 76952, the dynamical 
system mass is consistent with the theoretical system mass (Figure 
[8] top panel). The ages of the systems estimated from the position 
of the primary on the CMD are consistent with their position within 
the Local Interstellar Bubble (LIB). The minimum age of a star 
within the LIB, excluding thos e with relatively hi gh space motions 
(e.g. the ft Pic moving group - lOrtega et all2002h . has been shown 
to be 160 Myrs lAbll201ll) . 

The dynamical system masses of the three remaining systems 
are consistently lower than their theoretical system masses. One 



possible explanation for the apparently low dynamical masses is 
a non-Solar metallicity. Varying the metallicity has a significant 
effect on the system age estimate and the mass-magnitude rela- 
tions. As an example, a 2 M© star with super-Solar metallicity will 
be more luminous and have a redder V — K colour index than a 
similar-mass star of Solar metallicity. A super-Solar metallicity star 
will appear to be significantly older based on its position on a CMD 
if Solar metallicity models are used. To explore the effect of our as- 
sumption of Solar metallicity for the entire sample, the HIP 5300 
system (Figure[8] second panel) was studied at varying metallicity 
values. Using the Solar metallicity models, the dynamical system 
mass of this system is significantly lower than the theoretical sys- 
tem mass (Figure [9] top panel). Increasing the metallicity causes 
the star to appear both younger, and less massive (Figure[9] bottom 
panel). Only eight of the targets included within this study have 
metallicity measurement s, either from spectroscopic analysis (e.g. 
Erspa mer & North] |2003l) or derived from Stromgren photometry 
(e.g. lSong et al]|200ll) . demonstrating the need for further study in 
this area. 

6. 1.2 K- and M-type binaries - HIP 44127 

The detection of two hierarchical systems, with pairs of lower-mass 
companions in a wide orbit around an A-type primary, allows for a 
comparison of the theoretical models in the low-mass regime where 
the models differ in the treatment of the contraction phase of these 
objects (e.g. Figure[3]l. Our observation of the HIP 44127 system, 
shown in FigureQJJ] resolves three components to this hierarchical 
system, with an A-star primary (A) separated by ~ 4" from two 
fainter, gravitationally bound, companions (BC). 

The orbit presented in the Sixth Orbit Catalog of the BC pair 
around the A-type primary is in disagreement with our recent ob- 
servations of this system. Although the phase coverage is insuffi- 
cient for a robust orbital determination, the high proper motion of 
the primary (Aa = —441.1 mas yr _1 , AS = —215.2 mas yr _1 ) 
suggests that the BC pair is co-moving. In addition, radial veloc- 
ity variations detected within the spectra of the primary suggest the 
presence of a spectr oscop ic component to this system with a period 
of ~ 1 1 years |Abt|[l965). The separation of this component was 
anticipated to be between 072 and Q','6, reaching maxim um sep- 
aration in the middle of 2007 jDocobo & Andraddl2006T) . We do 
not resolve any companion in any epoch of our observations which 
are consistent with these predictions, placing an upper limit to the 
separation of 0'.'08, O'.'lO, and O'.'lO in 2008, 2010, and 2011 re- 
spectively. The magnitude difference between the primary and the 
suggested companion, estimated to be Am = 1.2, is well within 
our detection limits at the expected separation jDocobo & An drade 
120061) . 

A refined orbit of HIP 44127 BC is presented in Tableland 
Figure [6] with an estimated dynamical system mass of 0.68±0.04 
Mq . The small magnitude difference between the two components 
(AK ~ 0) suggests a mass ratio close to unity. Assuming indi- 
vidual masses of ~ 0. 34 Mp the components w ould be of early- 
to mid-M spectral type (Baraffe & Chabrier 1996), a region of par- 
ticular disa greement between theoretical models. The lBaraffe et al.l 
( 1998) and Siess et al.l feOOd) models both predict that for a binary 
consisting of two stars of magnitudes equal to the measured magni- 
tudes of the two components (Mk = 5.87), the system mass will 
increase from 0.3 to 0.9 M between 10 and 100 Myrs (FigurefTOt, 
as these models take into account the phase of contraction onto the 
Main Sequence for l ower-m ass stars. No such change is predicted 
by the lMarigo et al. I J2008h models, with the system mass remain- 
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Figure 10. Top row - (left): An AO observation of HIP 44127 triple system. The orbit of the lower-mass pair (BC) taken from Table[8]is plotted as a solid 
line for reference. The bright point source to the lower-right of the primary PSF is a known filter ghost. The axis labels are given in arcseconds, with each 
arcseco nd equal to approxim ately 14.5 AU. (centre): A c olour-magnitude diagr am highlighting the p osition of the primary with respect to the theoretical 
models: lMarigo et all 120081) - solid line. lSiess etaD l200fj|) - dotted line, while the Baraffe et al. 1 1998) models do not extend to high enough masses. The age 
of the system was estimated using the isochrone which best fit the position of the primary, (right): A colour-magnitude diagram highlighting the position of 
the two lower-mass components (B - blue, C - red). The line styles are the same as the centre panel Bottom row - The theoretical system mass plotted as a 
function of system age for each of the three models. The vertical dotted line within each panel indicates the age of the syste m estimated from using each of 
the models, while the horizontal dot-dashed line indicates the dynamical system mass estimated from the orbital fit. As the Bar affe et alj 119981) models do not 
extend to the A-type star mass range, an age estimate is not possible. 



ing unchanged between 10 Myrs and 1 Gyr (Figure [Tot. The theo- 
retical system masses are only cons i stent w ith t he dynamic a l mass 
at ~ 30 Myrs in the iBaraffe etaD dl998l) and ISiess et ail j2OO0t) 
model s, and are not consistent with any age in the iMarigo et al.l 
(2008) models. The metallicity o f this system has been measured to 
be almost Solar dWu et al' 2010), removing the degeneracy which 
exists between metallicity and the mass and age of the system (e.g. 
§6.2.6). 

The position of the primary on the CMD suggests a relatively 
young age of the system between 50 and 250 Myr. The V — K 
colours of the B C pair also suggest a y oung age, between 4 and 
1 00 Myr using the lSiess et aD fcOOCh and lBaraffe etail d 19981) mod- 
els. The inferred young age of the system is not inconsiste nt with 
the minimum age of stars found within the LIB (Abt 201 ll) . given 
the relatively high UVW space velocity of the HIP 44127 system 
dPalous|[l983l) . The detection of X-ray emission from this system 
is also of interest. Previous studies have shown that A7 stars such 
as the primary should not emit X-rays, and that any detection of 
X-rays from the positi on of the star can be indicative of a lower- 
mass companion (e.g. |Pe Rosa et ail l201ll) . The ROSAT source 



J08591 3.0+480227 is co incident with the optical position of HIP 
44127 i lVoges et all 19991) . 

This hierarchical triple system warrants further study, primar- 
ily in order to refine the orbital fit as the BC pair approaches apas- 
tron passage in 2018. This system also makes an ideal candidate 
for future spectroscopic observations to search for the narrow spec- 
tral lines of the two faint companions. With a double-lined spec- 
troscopic orbit fit, model-independent masses can be calculated for 
the individual components. If the young age suggested by both the 
position on the CMD and the system mass estimated from the mass- 
magnitude relations is correct, the BC pair would be an ideal cali- 
brator for the theoretical models in the young, low-mass regime. 

6. 1.3 K- and M-type binaries - HIP 82321 

The hierarchical triple system HIP 82321 was resolved within a sin- 
gle epoch of our AO observations. The A-type primary (A) is ~ 2" 
from a binary pair of two lower-mass companions (BC) in a wide 
orbit. The significant proper motion of the primary (Aq = 22.8 
mas yr~ , AS — —51.4 mas yr _1 ), and near constant separation 
of the three stars, suggest that the lower-mass pair is co-moving 
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Figure 11. The theoretical system mass plotted as a function of system age for each of the three models for the HIP 82321 BC system. The vertical dotted line 
within each panel indicates the age of the system estimated from using each of the models, while the horizontal dot-dashed line indicates the dynamical system 
mass estimated from the orbital fit. The clear discrepancy between the dynamical and theoretical system masses can be explained by the lack of complete 
coverage of the orbit presented in Figure|6] Further high angular resolution observations of this system will allow for an improved determination of the orbital 
parameters. 



with the primary. Spectroscopic measurements of the lower-mass 
components of this system are possible given the small magnitude 
difference of the pair with respect to the primary ( AK ~ 2.5); such 
spectra would allow for a determination of the individual masses 
independent of the distance to the system. The Hipparcos paral- 
lax and 2MASS A"-band magnitude measurements of this system 
allow for a tight constraint of the system age to between 300 and 
400Myrs, based on the position of the primary on the CMD. The 
primary is also a possible member of the Ursa Major moving group, 
an ass ociation of Solar metallicit y stars with an age between 30 
Myrs iSoderblom & Mavorll993T) and 500 Myrs fang et al.l2003h . 

A refined orbit of the BC pair is presented in Table[8] with an 
estimated dynamical system mass of 1.32±0.05 M Q . The theoret- 
ical system mass, assuming the distance to each component is the 
same as the primary, is significantly higher at ~ 2.0 Mq (Figure 
II U . The sparse coverage of the orbit (Figure [6]( suggests that the 
orbit fit could be poor, resulting in an incorrect dynamical system 
mass. The orbit fit would be significantly improved with subsequent 
observations. An alternative scenario is that the pair are a back- 
ground binary with a proper motion similar to the primary, which 
would bring the dynamical and photometrically derived masses into 
agreement. 

6.2 Higher-order multiplicity 

6.2.1 Known spectroscopic binaries 

The presence of an unresolved spectroscopic binary can have a sig- 
nificant impact on the magnitude assigned to each component of 
the multiple system, increasing the estimated component magni- 
tude by as much as 0.75 mag for an unresolved equal-mass spec- 
troscopic binary. For each target within this study, the literature was 
searched for references to additional components resolved through 
spectroscopic or interferometric observation which would influence 
any comparison made between the dynamical system mass and the 
theoretical system mass (e.g. Figure [T3l>. Of the targets within the 
orbit subsample, three are known to have additional spectroscopic 
components (HIPs 28614, 36850 , 80628). T he spectroscopic com- 
ponent to HIP 9480 resolved bv lAbtl (l965) is resolved within our 
adaptive optics observations. Similarly, for the monitoring subsam- 
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Figure 12. A schematic of the 7 hierarchical systems within both subsam- 
ples, including only components within a projected separation of 100AU 
to the A-type primary. The schematic diagram does not include the sus- 
pected spectroscopic binaries described in §6.2.2. Those systems which 
were newly resolved as a part of the VAST survey are highlighted in blue. 
Each filled circle represents an individual component, with the period of 
the binary pair denoted by its vertical position. Our high resolution obser- 
vations are only sensitive to systems with orbital periods greater than ~ 10 
years, denoted by the dashed horizontal line. For wide separation systems 
with indeterminate orbital periods, such as HIP 128 AB-C, the orbital pe- 
riod has been arbitrarily set at 1000 years, denoted by the dotted horizontal 
line. 



pie, two are known spectroscopic binaries (HIPs 128, 1 1661 1). The 
greater frequency of spectroscopic binaries within the orbit sample, 
relative to the monitoring sample, can be explained by the narrow 
spectral lines of the former sample due to their relatively low radial 
velocities, (vsini) = 71 km s ~ compared with (v sin i) = 115 
km s _1 for the latter sample dAbt & Morrelllll995h . The magni- 
tudes listed for the resolved components of these systems (Tabled 
are the blended magnitudes of the listed spectroscopic components. 
In addition to these known spectroscopic binaries, three members 
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Figure 13. The six systems for which spectroscopic components are either known to be present through previous spectroscopic observations (HIPs 28614, 
36850, 80628), or thought to exist due to a significant discrepanc y between the theoretical s ystem mass and the dynamical system mass. For clarity only 
one model has been used to estimate the theoretical system mass iLeie une & Schaerer 2001). The horizontal displacement between the intersection of the 
dynamical system mass (dot-dashed line) and the age estimate of the system (dotted line), and the theoretical system mass curve (solid shaded region), gives 
an order of magnitude estimate of the mass of the unresolved spectroscopic components. 



of the orbit subsample are resolved as hierarchical triples within 
our high resolution observations (HIPs 11569, 44128, 82321). A 
schematic representation of the higher-order multiplicity systems 
is given in FigurefPH 

For the three systems with known spectroscopic components 
within the orbit sample, the dynamical system mass includes the 
mass of each component, regardless of whether it is resolved within 
our data. This causes a significant discrepancy when the dynamical 
system mass is compared with the theoretical system mass (e.g. 
Figure I7J, with the dynamical system mass being systematically 
higher. The discrepancy between the two values cannot be directly 
converted into a mass for the unresolved components however, as 
a blended magnitude would have been used when determining the 
mass from the mass-magnitude relations. 



6.2.2 Evidence of spectroscopic components 

Systems with significantly higher dynamical masses than theoreti- 
cal masses obtained from mass-magnitude relations are strong can- 
didates for multiple systems with unresolved components which 
have not been detected in spectroscopic observations. We find two 
such systems with the signature of a possible unresolved compo- 
nent: HIP 17954 and HIP 93506 (Figure[[3j, with masses of the or- 
der of 0.5 and 1 .5 Mq . Sensitive spectroscopic observations of both 



systems may lead to the detection of the spectral lines from an unre- 
solved lower-mass component. A similar phenomenon is observed 
for the HIP 47479 system, although the number of measurements 
used to determine the orbit is particularly low, making the dynam- 
ical mass less certain. The narrow spectral lines of HI P 47479, im- 
plied by the low measured stellar rotational velocity jRover et al.l 
2007), make it an ideal candidate for spectroscopic follow up. Pre- 
vious spectroscopic observations of thi s system reve al v sin i vari- 
ations with a magnitude of 40 km s - 1 dMoordl 1932b . 



From the total sample of 26 systems, and only considering 
stellar companions within 100 AU, a lower-limit on the higher- 
order multiplicity of A-type stars can be estimated. Assuming the 
suspected unresolved companions described earlier in this section 
are true, there are five double, six triple, and two quadruple systems 
within the orbit subsample, corresponding to frequencies of 39, 46, 
and 15 per cent, respectively. This lower-limit shows an enhance- 
ment on the higher-order multiplicity of A-type stars when com- 
pared with Solar-type primaries (74 per cent dou ble, 20 per cent 
triple , and 6 per cent quadruple or higher-order - iRaghavan et al.l 
2010), and is more consistent with the fraction reported for more 
massive O-type primaries (46 per ce nt double, and 54 per cent triple 
or higher-order - Mason et ah 20091). 
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Figure 14. Observations of twelve of the thirteen systems with projected separations < 100AU suggested as future orbital monitoring targets. The companion 
in each image is highlighted, with a logarithmic image scale between 1 (white) and 10 -6 (black) relative to the peak intensity of the primary. 
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Figure 14. f continued): An observation of one of the thirteen systems with projected separations < 100AU suggested as future orbital monitoring targets. The 
companion in each image is highlighted, with a logarithmic image scale between 1 (white) and 10~ 6 (black) relative to the peak intensity of the primary. 



6.3 Continued monitoring targets 

6.3.1 Newly-resolved binaries 

We have identified 13 binary systems with projected separations 
ranging between 13 AU and 96 AU which would make ideal can- 
didates for future orbital monitoring projects (Table [5}. The 100 
AU projected separation cut-off was applied to select only systems 
for which orbital motion could be detected with several years of 
observations. The binaries resolved within the monitoring subsam- 
ple typically have lower-mass ratios than for the orbit subsample, 
a demonstration of the effectiveness of AO observations at detect- 
ing high-contrast binaries. Based on their position on the CMD, 
two members of the monitoring subsample (HIP 5310, HIP 18217) 
appear to lie on the Zero-Age Main Sequence, and the measured 
magnitude difference between primary and secondary would corre- 
spond to a late K or early M-type companion in each case. These 
companions are of particular interest as they will allow for further 
tests of the theoretical models within the young, low-mass regime. 
A gallery of the observations obtained of each target within the 
monitoring subsample is shown in FigurJPfl 



6.3.2 HIP 77660 

An advantage of monitoring binary systems for orbital motion with 
AO instrumentation over interferometric techniques is the elimina- 
tion of the quadrant ambiguity. In some cases, the output of the 
image processing of speckle dat a results in a 180° am biguity on 
the position angle measurement ( Bagn uoloet ai1ll992t) . This may 
lead to a scenario where orbital motion is thought to exist for a bi- 
nary pair, when in fact this motion is the product of such a quadrant 
uncertainty combined with the true linear motion of the compan- 
ion. The ambiguity can be resolved by observing the system using 
AO imaging, where no reconstruction is required to obtain the final 
science image. 

Based on our AO observations of the binary system HIP 
77660, it appears that such a quadrant ambiguity has occurred, 
making linear motion appear as orbital motion, as shown in Fig- 
ure [T5] The significant proper motion of the system (Aa = —91 
mas/yr, AS — — 28 mas/yr), as measured by Hipparcos, is incon- 
sistent with a stationary background object. While future measure- 
ments of this system will be able to resolve the presence of or- 
bital motion after a sufficient time baseline, spectroscopy or multi- 



colour photometry will allow for a rapid characterisation of the 
properties of both components. 



7 SUMMARY 

We have presented high resolution observations of 26 nearby mul- 
tiple systems with A-type primaries with projected separations < 
100 AU, 11 of which are binaries newly resolved as a part of the 
VAST survey. For those systems with sufficient orbital motion cov- 
erage, refined orbital parameters were calculated and the estimated 
dynamical system mass was compared to masses derived from the- 
oretical models. Due to their rapid evolution across the CMD, re- 
moving the significant age degeneracy for lower-mass solar-type 
stars, binaries with A-type components are ideal targets with which 
to test theoretical models. Four such systems were investigated, 
with one system having consistent dynamical and theoretical sys- 
tem mass estimates. Of the remaining three systems, each had a 
dynamical mass significantly lower than that predicted from the 
models. While this discrepancy may be indicative of a true diver- 
gence between the models and observations, the lack of metallicity 
measurements for these systems provide another explanation. Fu- 
ture orbital monitoring observations of A-star binary systems will 
provide further refinement to the orbital parameters and, combined 
with refinement of the magnitude, metallicity, and parallax mea- 
surements, will improve the analysis performed within this study. 

Observations of two hierarchical systems, consisting of an A- 
type primary and a tight low-mass binary pair in a wide orbit al- 
lowed us to extend this analysis to the lower mass regime. The pri- 
mary of the triple HIP 44127 suggests a young system age (<100 
Myr), making the M-dwarf pair interesting for comparison with 
evolutionary models. We have also shown that a dynamical system 
mass significantly higher than the theoretical system mass is sug- 
gestive of an unresolved spectroscopic component within the sys- 
tem. Demonstrated on several known spectroscopic binaries, sys- 
tems which exhibit this discrepancy are ideal candidates for fu- 
ture spectroscopic or interferometric observations in an attempt to 
detect these hypothetical companions. Interferometric observations 
may be required to resolve these purported companions, as the ro- 
tationally broadened spectr al lines of the r apidly rotating members 
of the monitoring sample jAbt & Morrell|[l995l) may preclude the 
detection of additional components using spectroscopy. Including 
the three systems with evidence of an unresolved close companion, 
a lower-limit on the higher-order multiplicity can be estimated from 
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Figure 15. (top panel): The historical measureme nt of the syst e m HIP 
77660, plotted with the orbital fit estimated by iDocobo et a?] J2010h . 
This system demonstrates the quadrant ambiguity limitation of the 
speckle interferometry technique. The binary companion was measured by 
iTokovimn etal] feOlfJ) using speckle interferometry to be at a position an- 
gle of 248.4° to 246.8° between mid-2007 and mid-2009 (blue filled cir- 
cles, top right). Our AO images, which do not suffer from the same quadrant 
ambiguity, shows the binary at a position angle of 69.7° to 67.3° between 
mid-2006 and mid-2011 (red filled stars, left), (bottom panel): The relative 
position of the secondary, with 180° added to the four speckle interferome- 
try measurements which were inconsistent with the observations presented 
within this study. The solid line connects the first and last observations, 
showing a deviation away from linear motion. Further observations are re- 
quired to determine if this deviation is indicative of the true orbital path of 
the secondary about the primary. 



the 13 systems within the orbit subsample as 39 per cent double, 46 
per cent triple, and 15 per cent quadruple. The frequency of single 
A-type stars will be explored in an upcoming publication within 
the VAST survey paper series. The remaining systems for which 
an orbit could not be determined are candidates for orbital moni- 
toring projects with ground-based high resolution observations. A 
number of these systems are of particular interest, based on age es- 
timates derived from the position of the primary on the CMD, and 
the magnitude difference between the two components. 
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